Structural properties of dense hard spheres near random close packing 
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The structural properties of dense random packings of identical hard spheres (HS) are investigated 
by using large scale simulations. The bond order parameter method is used to obtain detailed 
information on the local orientational order of the HS system at different packing fractions (j>, in 
the range between <j> « 0.53 and <f> « 0.72. The probability distribution functions of the spheres 
versus different qi wi are used to construct new measures, characterizing the crystallization of the 
system. Our results show that behind the Bernal limit the crystalline clusters transform into the 
global three-dimensional crystalline structure, which, upon further densification, transforms into 
alternating planar layers formed by different lattice types. 

The model of hard spheres (HS) is of fundamental importance in condensed matter physics and material science 
since it successfully reproduces the essential structural properties of liquids, crystals, glasses, colloidal suspensions 
and granular media [J 0. Structural changes of the HS system have been observed when a packing is densified 
above the density of random close packing (RCP) <j) c ~ 0.64 (see, e.g. [!, H|). A new order parameter, based on 
the cumulative properties of spheres distribution over the rotational invariant w§, was proposed in Q to identify this 
structural changes. Random arrangement is transferred to partial crystallization across RCP, and two lattice types 
(fee and hep) are observed in the crystalline clusters. However, more detailed information has not been obtained due to 
insufficient system size. In this study, we aim to test previous findings, and to further explore the structural properties 
of the transition at RCP in more details in comparison with previous studies [3, 0], based on larger scale simulations. 
In particular, we take a large set of packings composed of N = 64 x 10 3 identical spheres with periodic boundary 
conditions, generated using modified Lubachevsky-Stillinger (LS) algorithm Q. The packing fraction depends on 
the compression rate as slower compressions result in denser packings. The packings obtained from LS algorithm 
are then used as initial configurations to generate mechanically stable jammed packings of particles interacting with 
Hertz normal forces. When (f> < 0.64, Mindlin frictional forces are added to stabilize the packings. More details of the 
numerical algorithm are described in Refs. [4[. The corresponding packing fractions vary in the range cj> ~ 0.53 — 0.72. 
The range includes the random close packing state at C ~ 0.64 (Bernal limit) [f|. 

To define the local structural properties of the system we use the bond order parameter method (3 , which has been 
widely used in the context of condensed matter physics HS systems [lOHlql. complex plasmas [16-21|, colloidal 

suspensions (2^, H3i granular media 24 1 etc. In this method the rotational invariants (RIs) of rank I of both second 
qi (i) and third wi (i) order are calculated for each sphere i in the system from the vectors (bonds) connecting its center 
with the centers of its N nn (i) nearest neighboring spheres: 
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centers of spheres i and ?'. In Eq.© are the Wigner 37-symbols, and the summation in the latter 
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expression is performed over all the indexes wii = —I, I satisfying the condition mi+m,2 + m,3 = 0. In 3D the densest 

possible packings of identical hard spheres are known to be face centered cubic (fee) and hexagonal close-packed (hep) 

lattices (4>f cc — 4>hcp = v27r/6 ~ 0.74). Presence of spheres having icosahedral (ico) symmetry is also possible in the 

dense HS system. To detect these structures we calculate the rotational invariants qi, qe, u>4 and we for each sphere 

using the fixed number of nearest neighbors N nn = 12 (25[. A sphere whose coordinates in the 4-dimensional space 

(<j4, qe,W4,we) are sufficiently close to those of the ideal fee (hep, ico) lattice is counted as fee-like (hep-like, ico-like) 

sphere. By calculating the bond order parameters it is easy to identify also disordered (liquid-like) phase (for instance, 



2 




FIG. 1: (Color online) Probability distribution functions (PDFs) of the rotational invariant qi (q§ (a) and q$ (c)) and Wi 
(wq (b) and W4 (d)) for different packing fraction <j>. Values of the bond order parameters for ideal lattice types are also 
indicated. Cumulative distributions (normalized to unity) are shown by dashed line of the same color; the distributions are 
much less noisy than PDFs. Insets (from top to bottom) show HS distributions over the cubic box at different packing fraction 
cj> ~ 0.65, 0.66, 0.68, respectively. The spheres are color-coded by q§ value to visualize disordered liquid-like (blue), hep- like 
(green) and fee- like (red) phases. Number of particles in the box is N — 64 x 10 3 . 
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FIG. 2: (Color online) Order parameter Wg h (red line and black squares) and the mean value of the rotational invariant (q&) 
(green line and blue squares) versus packing fraction <j>. Lines and circles correspond to HS system with N = 10 4 spheres; 
squares correspond to N = 64 x 10 3 spheres. Good agreements between different system sizes are clearly seen. Explosive-like 
growth of the parameter Wg h ((f)) at (j> > (f> c ~ 0.64 is a signature of the appearance of crystalline order (hep-like and fcc-like 
spheres) in the system. On the other hand, increase in (qe) near <j> c is monotonous and relatively slow. The possibility for the 
dense HS systems to have quite different compositions (with different relative density of fee and hep crystalline configurations) 
is the main reason behind the pronounceable oscillations in (qa) observed for <f> > 0.68. Inset shows number of ico-like spheres 
versus <f>: ico-like spheres exist only in the narrow range of parameter <f) near RCP transition. 



such spheres have mean bond order parameter g 6 iq ~ N n n^ 2 ~ 0.29 <C q f £ c / hc P/ lco ^ b v varying number of nearest 
neighbors N nn and rank I of bond order parameter it is possible to identify any lattice type (including quasicrystalline 
particles and distorted hcp/fcc/ico modifications) existing in the system of particles (e.g, by using N nn — 8 and 14 it 
is easily to identify the first and the second shells of the body centered cubic (bec) lattices, etc). We note, that the 
bond order parameters wi oc qf, so, in general, the parameters wi are much more sensitive to the local orientational 
order in comparison with qi . 

Figure [T] shows the probability distribution function (PDF) of the rotational invariants of second q^ q§ and third 
u>4, u>6 orders, for different packing fraction <p covering both liquid-like and crystalline states of HS. Densification of 
the HS, as clearly seen in the Figure [TJ results in appearance of the crystalline particles in the vicinity of the Bernal 
limit (RCP). Cumulative distributions (shown by dashed lines) are also plotted; those can be used to define the 
number density of hcp/fcc/ico- like spheres in the system (and disordered solid- like spheres as well). Such cumulative 
distributions being much more smoothed than corresponding PDFs are proposed to create the order parameters 
characterizing the phase state of the HS system. 

It has been shown recently (2(3 | that the cumulative distribution function of the form 

W 6 (x)= f n(w 6 )dw 6 (3) 

is an extremely sensitive measure of the structural order in the system (in Ref. [2(| the function Wq(x) was introduced 
for the first time to quantify the solid-liquid transition in confined 3D Yukawa system). Here n(we) is the PDF of 
spheres over the rotational invariant we normalized to unity. The relevant order parameter is the position of the half- 
height of the cumulative distribution u>g h , such that W^Wq 11 ) = 1/2. Here, the order parameter is used to identify 
the onset of crystallization in the randomly packed HS system. 

Figure [2] shows the dependence of the order parameter u>g h and the mean value of [defined as (q 6 ) = 
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N 

(1/iV) on the packing fraction <fi. Lines and circles corresponds to HS system with N = 10 4 spheres [3j; 

i=l 

squares corresponds to the case with TV — 64 x 10 3 spheres. The very good agreements between different system sizes 
are clearly seen — it suggests that HS consisting of N — 10 4 spheres well reproduces the global structural properties 
of the HS system. Explosive-like growth of the order parameter Wg h at (f> > </> c reflects the appearance of crystalline 
order (hep- like and fcc-like spheres). This proves the usefulness of this structural measure to locate the onset of crys- 
tallization, as demonstrated here in the context of randomly packed HS systems. The dependence (qq) on <fi shown 
in Fig. [5] also demonstrates interesting properties. Relatively slow growth of (q 6 ) for 4> ^ 4>c is associated with the 
appearance of single hep-like spheres [3|- For 4> > 4>c the parameter (q 6 ) grows more rapidly. Further densification of 
the system results in strong oscillations of (qs), which start at (f> ~ 0.68 and suggest another transition Geometry 
gives a simple explanation of such an oscillatory behavior: The point is that the hep and fee crystal structures are 
characterized by quite different values of the parameter q 6 (for ideal hep and fee lattices <7g Cp ss 0.48 and q f g c « 0.57). 
Dense HS packings can be apparently realized (for a given 4>) with significantly different relative quantity of hep and 
fee phases. This explains the oscillations. Inset in Fig. [2] shows another important measure characterizing structure 
of HS near RCP: number density N- lco of spheres having icosahedral symmetry (ico-like or fivefold spheres) versus 
packing fraction <p. It shows that ico-like spheres can exist only in a narrow range of parameter <f>: at low <p number 
of ico-like spheres is exponentially small; densification of HS above <f> ~ 0.66 results in nearly disappearance of the 
ico-like spheres. This suggests that the abundance of ico-like spheres can be used as a measure, characterizing HS 
crystallization. 

Cluster analysis of the arrangements of hep-like and fcc-like spheres reveals another interesting structural property 
of dense HS packings. Figure [3] shows how the distribution and the shape of clusters composed of hep and fcc-like 
spheres vary upon the densification of the system. For relatively low volume fractions <fi < 0.65 (Fig.[3fa, b)), fcc-like 
clusters consist of only few spheres. Increase of (j) results in the appearance of several big (containing few tens of 
spheres) both hep and fcc-like clusters. A typical cluster has a complicated 3D shape (Fig. [3jc, d) for <j> rj 0.66). 
Eventually, nuclei of fcc-like and hep-like clusters transform into a global crystalline structure occupying most of 
the system volume. Further increase in <fi leads to a structural change: a transition from 3D clusters to 2D layers 
composed of fcc-like and hep- like spheres occurs which is in full agreement with results of 

The relative numbers of spheres in hep-like (nhep) and fcc-like (nf cc ) states are plotted in Figure [4] as functions of 
<f>. For 4> < 4> c the relatively slow increase of n to t = ^hep + nfcc with (f> is associated with emergence of single hep-like 
spheres, which are nearly randomly distribute over the system volume. At ~ cj) c the slopes of n t ot(</>) and rif cc ((j>) 
are drastically increase. The number of fcc-like spheres reaches that of the hep-like spheres (nf cc ss nhep) &t 4> — 0.66, 
where transition from hep-dominated to fee-dominated packing occurs. Additionally, the abundance of spheres having 
number of contacts N c = 12 is plotted in Figure [4j making this measure to be one more very sensitive indicator of 
the crystallization of a dense HS. Inset in Figure 0] shows the relative cumulative spectrum of hep (fcc)-like (red and 
green line, respectively ) clusters as a function of the cluster mass at different (j) values. The hep/fee spectrum near 
the RCP ((f) ~ 0.65) reveals already mentioned features: only small clusters consisting of few spheres are present in 
the HS system. As the density increases, larger clusters appear. This observation is consistent with the classical 
nucleation process 26] ; the denser packings are more "supercooled" before falling out of equilibrium, therefore they 
have lower free energy barriers to form large clusters. 

The above results illustrate a sharp first-order like transition [ij at 4> c between disordered (frictional) packings 
and partially ordered packings. Amorphous frictionless packings only exist at a single density <j) c : below C , the 
packings are unstable unless they are frictional (and the frictional packings are disordered); above (j> c , the packings 
are partially ordered. In contrast, multiple amorphous frictionless states have been found in other simulations (27l.l28j. 
which supports the mean-field theoretical prediction that amorphous jammed packings should exist in the interval (so- 
called J-line) <fi £ [4>th, </>gcp] The key point here is whether crystallization is present or not. If the nucleation rate 



is slow enough, such as in polydisperse |27j or large dimensional systems [28], then crystallization is suppressed and the 
J-line is well observed. In the present case we have employed monodisperse sphere packings in three dimensions. In 
such systems, the nucleation rate is relatively fast, which prevents any slow compression to produce dense amorphous 
packings. As a result, we find any frictionless packings above <p c are partially crystallized as revealed by measures 
of order parameters. Combined with the mean-field prediction, our results suggest two possible scenarios in 3D 
monodisperse packings: (i) <p c < 0th, therefore amorphous packings are hidden by partial crystallization; (ii) t h ^ </>c, 
however the difference between t h and <f> c is small and difficult to detect based on present numerical accuracy. Further 
test of the subtle differences between these two scenarios may require even larger-scale simulations. 

To conclude, we have demonstrated that the order parameter, based on the cumulative distribution of spheres over 
the values of the rotational invariant wq is a very sensitive measure to describe the onset of crystallization in the 
system of randomly packed hard spheres. The proposed order parameter is also especially convenient indicator of the 
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appearance of hep , fee and ico-like crystalline particles in the system. It has been shown, that abundance of icosa- 
hedral particles decrease drastically after RCP limit, making the measure to be one more indicator of the HS system 
crystallization. We used this to investigate how the relative distribution between hep and fee arrangements varies 
with increasing the sphere packing fraction. Finally, we observed that at <j) — 0-68 the structural transition occurs: 
The essentially three dimensional shape of the global crystalline aggregate changes to planar (layered) structure. 

We thank F. Zamponi and P. Charbonneau for helpful discussions. This study was supported by NSF-CMMT and 
DOE Geosciences Division. B. K. was partly supported by the Russian Foundation for Basic Research, Projects No. 
13-02-00913 and 13-02-01099. 
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FIG. 3: (Color online) Distribution of hep and fee clusters over the space at different <f> values: <j> ~ 0.65 (a, b), (j> ~ 0.66 (c, d) 
and 4> — 0-68 (e, f). Particles are color-coded by the mass of the cluster (in the unit of the mass of a single sphere). At low cj> 
single hep-like spheres dominate (a, b). The densification leads to the formation of large 3D clusters of both hep and fee- like 
spheres. Eventually, clusters of fcc-like and hep-like spheres transform into a global crystalline structure. Further densification 
results in the structural transition, when essentially 3D shape of this global crystalline structure changes to planar 2D layers 
formed by fcc-like (shown in the figure) and hep-like spheres. 
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FIG. 4: (Color online) Abundance of hep-like (red) and fcc-like (green) spheres versus the packing fraction <j> as seen from 
detailed simulations of HS systems. For moderate packing fractions (<f> < <f> c ), the increase of crystalline particle with (f> is 
due to the appearance of a single hep-like spheres. Results of large scale simulations with N = 64 x 10 3 spheres are shown 
by squares (orange and green colors correspond to hep and fcc-like spheres). The dependencies show pretty well agreement 
between the cited HS systems. Additionally, the abundance of spheres having number of contacts iV c = 12 is plotted by blue 
square, making this measure to be one more sensitive indicator of the crystallization of dense HS systems. Inset shows the 
relative cumulative spectrum of hep (fcc)-like (red and green line, respectively ) clusters versus its mass at different (f> values: 
<fi ~ 0.68 (solid, triangles), <j> ~ 0.66 (dashed, gradients) and cf> ~ 0.65 (solid, circles). 



